home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / dlags2.z / dlags2
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. DDDDLLLLAAAAGGGGSSSS2222((((3333FFFF))))                                                          DDDDLLLLAAAAGGGGSSSS2222((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DLAGS2 - compute 2-by-2 orthogonal matrices U, V and Q, such that if (
  10.      UPPER ) then   U'*A*Q = U'*( A1 A2 )*Q = ( x 0 )  ( 0 A3 ) ( x x ) and
  11.      V'*B*Q = V'*( B1 B2 )*Q = ( x 0 )  ( 0 B3 ) ( x x )  or if ( .NOT.UPPER )
  12.      then   U'*A*Q = U'*( A1 0 )*Q = ( x x )  ( A2 A3 ) ( 0 x ) and  V'*B*Q =
  13.      V'*( B1 0 )*Q = ( x x )  ( B2 B3 ) ( 0 x )  The rows of the transformed A
  14.      and B are parallel, where   U = ( CSU SNU ), V = ( CSV SNV ), Q = ( CSQ
  15.      SNQ )  ( -SNU CSU ) ( -SNV CSV ) ( -SNQ CSQ )  Z' denotes the transpose
  16.      of Z
  17.  
  18. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  19.      SUBROUTINE DLAGS2( UPPER, A1, A2, A3, B1, B2, B3, CSU, SNU, CSV, SNV,
  20.                         CSQ, SNQ )
  21.  
  22.          LOGICAL        UPPER
  23.  
  24.          DOUBLE         PRECISION A1, A2, A3, B1, B2, B3, CSQ, CSU, CSV, SNQ,
  25.                         SNU, SNV
  26.  
  27. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  28.      DLAGS2 computes 2-by-2 orthogonal matrices U, V and Q, such that if (
  29.      UPPER ) then
  30.  
  31.  
  32.  
  33. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  34.      UPPER   (input) LOGICAL
  35.              = .TRUE.: the input matrices A and B are upper triangular.
  36.              = .FALSE.: the input matrices A and B are lower triangular.
  37.  
  38.      A1      (input) DOUBLE PRECISION
  39.              A2      (input) DOUBLE PRECISION A3      (input) DOUBLE PRECISION
  40.              On entry, A1, A2 and A3 are elements of the input 2-by-2 upper
  41.              (lower) triangular matrix A.
  42.  
  43.      B1      (input) DOUBLE PRECISION
  44.              B2      (input) DOUBLE PRECISION B3      (input) DOUBLE PRECISION
  45.              On entry, B1, B2 and B3 are elements of the input 2-by-2 upper
  46.              (lower) triangular matrix B.
  47.  
  48.      CSU     (output) DOUBLE PRECISION
  49.              SNU     (output) DOUBLE PRECISION The desired orthogonal matrix
  50.              U.
  51.  
  52.      CSV     (output) DOUBLE PRECISION
  53.              SNV     (output) DOUBLE PRECISION The desired orthogonal matrix
  54.              V.
  55.  
  56.      CSQ     (output) DOUBLE PRECISION
  57.              SNQ     (output) DOUBLE PRECISION The desired orthogonal matrix
  58.              Q.
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.